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ABSTRACT 

The fraction of quasar's and gamma-ray burst (GRB) afterglows that vary due to 
microlensing by the stellar populations of intervening elliptical/SO galaxies is computed 
by combining the joint distribution of effective microlensing convergence (k) and shear 
(7) with microlensing magnification patterns. Microlensing is common in multiply im- 
aged sources. We find that 1 in 3 multiply imaged quasars should vary by more than 0.5 
magnitudes per decade due to microlensing, while 10% of macrolensed GRB afterglows 
should show a departure of more than 0.5 magnitudes from their intrinsic light-curve 
during the first 30 days. However microlensing by stars is rare in general, with only 
1 source in ~ 500 varying by more than 0.5 magnitudes during the same periods. We 
find that most microlensing by stars will be observed in a regime where 7 > 0.1. Thus 
point-mass lenses do not provide an adequate description for most microlensing events. 
If dark matter halos contain a large fraction of mass in compact objects, the fraction 
of microlensed (by 0.5 magnitudes) images rises significantly to ~ 1 in 10 for quasars 
and ~1 in 5 for GRB afterglows. Comparison of variability between macrolensed and 
normal quasar images, and a moderate number of well sampled GRB afterglow light- 
curves should therefore discover or refute the existence of stellar mass compact objects 
in galaxy halos. While microlensing results in departures of the distribution of magnifi- 
cations from that of a smooth model, the effect on the macrolensing magnification bias 
for the discovery of lenses in quasar surveys is small. On the other hand, microlensing 
significantly broadens the distribution of macrolensed image flux ratios. 

Subject headings: gravitational lenses: microlensing - dark matter 



1. Introduction 

Cosmological microlensing of background quasars by stars in foreground galaxies was first 
discussed in the work of Chang & Refsdal (1979), Gott (1981) and Young (1981). Since then 
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considerable theoretical and observational progress has been made towards using microlensing phe- 
nomenon to obtain otherwise inaccessible information on the cosmological distribution of stellar 
and sub-stellar mass compact objects, and on the smallest scales of the central engines of quasars. 
In contrast, there has been relatively little work (Bartelmann & Schneider 1990; Koopmans &: 
Wambsganss 2000; Wyithe & Turner 2002, hereafter Paper I) extending calculations of the a-priori 
probability for cosmological microlensing beyond the seminal work of Press k, Gunn (1973). 

Microlensing was first identified in Q2237-I-0305 (Irwin et al. 1989; Corrigan et al. 1990). 
Early models presented by Wambsganss, Paczynski & Katz (1989) and Wambsganss, Paczynski 
&; Schneider (1990) explained the observed flux variation using a model containing stellar mass 
objects. Following this success, microlensing of cosmological sources (mostly quasars) has been used 
to discuss the nature of the intervening compact object populations (e.g. Press & Gunn 1973; Blaes 
& Webster 1992; Schneider 1993; Dalcanton et al. 1994; Lewis & Irwin 1995; Perna & Loeb 1997; 
Schmidt & Wambsganss 1998; Marani et al. 1999; Wyithe, Webster, Turner & Mortlock 2000; 
Wambsganss et al. 2000; Refsdal, StabeU, Pelt & Schild 2000; Schild 1999; Mao & Loeb 2001; 
Lewis & Ibata 2001; Lewis, Ibata & Wyithe 2000). Microlensing has also been used to constrain 
different physical emission regions of quasars. Variability of the optical continuum in Q2237-I-0305 
(Corrigan et al. 1989; Irwin et al. 1990; 0stensen et al. 1996; Wosniak et a. 2001a, b) has been used 
to show that the emission region is < lO^^cm (Wambsganss, Paczynski &: Schneider 1992; Ranch 
& Blandford 1992; Jaroszynski, Paczynski & Schneider 1992; Wyithe, Webster & Turner 2000). 
Similarly, data for Q0957+561 (Pelt, Schild, Refsdal & StabeU 1998) was used to find a source 
smaller than 5 x 10^"'' cm (Refsdal, StabeU, Pelt k. Schild 2000). On larger scales, Lewis et al. (1998) 
observed microlensing induced variation between the optical continuum and broad-line regions, and 
Wyithe, Agol & Fluke (2002) used mid-IR to V-band flux ratios to demonstrate that the mid-IR 
emission comes from a region > 10^ ''cm (based on the mid-IR data of Agol, Jones & Blaes 1999). 
While these results are specific to Q2237-I-0305 (particularly) and Q0957-I-561 both of which have 
long monitoring histories, monitoring programs for the determination of time-delays (e.g. Kundic 
et al. 1997; Schechter et al. 1997; Burud et al. 2000) in multiply imaged quasars offer a promising 
avenue for more microlensing observations. For example 0stensen et al. (1997) found evidence for 
low-level microlensing in monitoring of the Clovcr-lcaf quasar. 

More recently, additional novel microlensing observations have been published. These include 
evidence for microlensing of radio emission substructure in B1600-I-434 (Koopmans et al. 2000; 
Koopmans & de Bruyn 2000), evidence for microlensing induced polarization variability in II1413-I-1143 
(Chae et al. 2001), and evidence for microlensing of the Fe Ka line in MG J0414-I-0534 (Chartas 
et al. 2001). Several other potential claims of AGN microlensing have also been published (e.g. 
Lewis, Robb & Ibata 1999; Webb et al. 2000; Torres, Gustavo & Eiroa 2002). 

In addition to the above lists of observations and interpretations, many studies have pointed 
out the potential value of detailed microlensing observations (particularly of caustic crossings). 

The most exciting are schemes to determine the emission spectrum of the central engine as a 
function of radius (e.g. Grieger, Kayser &; Refsdal 1988; Grieger, Kayser & Schramm 1991; Agol & 
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Krolik 1999; Mineshige & Yonehara 1999) from multiband observations of caustic crossings. In other 

examples, Schneider & Wambsganss (1990) showed that the structure and kinematics of Broadhne 
regions could be constrained using observations of microlcnsed quasars, Belle &: Lewis (2000) have 
demonstrated that microlensing should result in the variability of polarized light, and Lewis, Ibata 
k. Wyithe (2000) and Lewis k, Ibata (2001) have shown that cosmological compact objects could 
be detected through surface brightness variations in giant cluster arcs and galaxies. 

Quasars are not the only cosmological sources subject to microlensing variability. Another 
candidate is the GRB afterglow (Loeb & Perna 1998). Indeed, Mao & Loeb (2001) find that all 
afterglows should be microlensed at a low level. Recently, a deviation from the generic double pow- 
erlaw decay was observed in monitoring of the afterglow light-curve of GRB000301C (Garnavich, 
Loeb &: Stanek 2000, based on data compiled by Sagar et al. 2000) and interpreted as due to 
microlensing. Detailed event inversions modeled after work on the interpretation of quasar caustic 
crossing events (e.g. Geiger et al. 1988) have since been performed (Gaudi, Granot k, Loeb 2001), 
and the lens hypothesis found to be plausible. Furthermore, Koopmans k Wambsganss (2000) have 
discussed the a-posteriori probability for microlensing of this source and found that while unlikely, 
the microlensing scenario cannot be ruled out. Their calculation went two steps beyond the formal- 
ism of Press k Gunn (1973). Firstly, they assumed microlenses to be clustered in galaxies which 
were modeled as SIS's of compact objects. Secondly, they computed microlensing statistics using 
magnification patterns rather than a point mass lens approximation. 

Microlensing of high redshift type la supernovae has been discussed by Metcalf k Silk (1999) 
and Wang (1999), and latter by Mortsell, Goodbar k Bergstrom (2001) and Minty, Heavens k 
Hawkins (2001) with emphasis on the utility of forthcoming survey samples. Type la supernovae 
are particularly interesting microlensing sources because of their status as standard candles. They 
are not expected to act as point sources for the typical duration of observations, and the resulting 
microlensing induced departures from the intrinsic light-curve will impact on the interpretation of 
the light-curve shape, and thus the inference of the intrinsic brightness. A study of microlensing 
and the light-curves of type la supernovae based on results described in this paper is currently in 
preparation and will be presented elsewhere. 

While microlensing promises to help with studies of cosmological sources and intervening com- 
pact object populations, it may also contribute to uncertainty in statistics of multiple imaging, both 
by changing the level of magnification bias in gravitational lens surveys and by causing variation in 

the flux ratios of multiply imaged quasars. Many authors have noted that optically measured flux 
ratios cannot be used as reliable constraints for galaxy lens models due to microlensing of one or 
more images (e.g. Schneider et al. 1988; Kochanek 1991; Witt, Mao k Schechter 1995; Mediavilla 
et al. 1998). Indeed, variation between optical and radio flux ratios (the radio should not be subject 
to microlensing) has been observed in several lensed quasars (e.g. Falco et al. 1996; Schechter k 
Moore 1993; Katz k Hewitt 1993). The effect of microlensing on magnification bias and flux ratios 
was discussed by Bartelmann k Schneider (1990). They suggested that the effect could be severe, 
particularly on the flux ratios. However, their model was simplistic. In particular they assumed 
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a point source, and analytic approximations for magnification patterns. On the other hand, with 
the exception of Koopmans & Wambsganss (2000) and Loeb & Perna (1997) (who used results 
of Bartelmann & Schneider 1990), this is to our knowledge the only previous paper to consider 
microlensing statistics based on the distribution of optical depth and shear for lines of sight to 
background cosmological sources. 

The growing base of different observations for an increasing number of systems motivates 
us to explore the a-priori probabilities for the observation of microlensing phenomena, as well as 
the circumstances under which that microlensing will be observed. The work is performed with 
views towards determining the types of microlensing observations that will be plausible in the 
future, and towards testing our understanding of the populations involved in much the same spirit 
as investigations of macrolensing probabilities (e.g. Turner, Ostrikcr &: Gott 1984; Turner 1990; 
Kochanek 1996). We begin in Sec. 2 by describing our model which is based on that described 
in Paper I, as well as the method for computing the probability distribution for the microlensing 
optical depth and shear. In Sec. 3 we discuss magnification distributions for lensed images including 
microlensing, and the cross-section for microlensing variability. Furthermore, we discuss the quasar 
microlensing variability statistics expected from two different classes of observational survey. We 
then discuss the effect of microlensing on the magnification bias for gravitational lens surveys and 
the distribution of image flux ratios for multiply imaged quasars. Finally in Sec. 4 we compute 
microlensing cross-sections for GRB afterglows. Macrolensing by galaxies results in two classes 
of images. Throughout the paper we refer to images of sources that are macrolensed as multiply 
imaged, and those that are not as singly imaged. We assume a flat cosmology having $7 = 0.3, 
A = 0.7 and Hq = 65 km sec^^ Mpc~^. Numerical tables of distributions presented in this paper 
will be made available upon request. 



2. The Joint Cross-Section for k and 7 

This section describes the calculation of the joint probability distribution for the microlens- 
ing optical depth k and shear 7 near images of randomly distributed cosmological sources. This 
probability distribution will be used in calculation of magnification distributions and microlensing 
statistics in subsequent sections. Paper I presented the distribution of k due to the normal stel- 
lar populations of galaxies for quasar images with random source positions. The model used is 
summarized below, however we refer the reader to Paper I for further details. 

A constant co-moving number-density of galaxies was assumed. The total mass-distributions of 
these galaxies were described by isothermal spheres having a central velocity dispersion ctdm = /c" 
(where a is the observed central velocity dispersion*^ and / ~ 1). The stellar populations embedded 



is the correction factor between the velocity dispersions of the luminous and dark matter (/ = y | obtained for 
the simplest dynamical models was introduced by Turner, Ostriker & Gott 1984). Kochanek (1996, and references 
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in these galaxies form the microlens population. In elliptical/SO galaxies, stars were distributed 

with dc-Vaucouleurs profiles (mass-to-light ratios were assumed constant in radius and red-shift), 
and in spiral galaxies with de-Vaucouleurs bulges and Kusmin (1956) discs. The distribution of 
velocity dispersions was described by the combination of a Schechter function d^/da and the 
Faber-Jackson (1976) relation. The characteristic velocity dispersion of an galaxy was taken to 
be (7^ = 220kmsec~^. 

In this work we follow a similar prescription. However, following the result of paper I that 
spiral galaxies contribute < 10% to the microlensing rate, we consider only elliptical/SO galaxies 
in our microlensing calculations. Furthermore, we allow the microlens mass-density to evolve with 
redshift in proportion to the cumulative star- formation history^ . A simple model which has constant 
star-formation from z = 10 till z = \ and then a rate proportional to (1 -|- z)'^ until the present 
day was used (e.g. Hogg 1999 and references therein; Nagamine, Cen &: Ostrikcr 2000). The 
mass-to-light ratios of the galaxies were then normalized so that the elliptical plus spiral galaxy 
populations contain the cosmological density in stars at redshift zero [fi^i, = 0.005 (Fukugita, Hogan 
& Peebles 1998)]. The model of a given elliptical galaxy with redshift Zg and central velocity 
dispersion a has an overall convergence at radius ^ of 

KSIS = ^ ;^7, (1) 



-'crit 



2Ge 



and an average local convergence in stars of 



-'crit 

where 



is the critical density for lensing, Rq and Sq (the characteristic radius and central density) are 
functions of a (see Paper I, after Djorgovski & Davis 1987), and T is the mass-to-light ratio of the 
stellar population. Assuming the dark halo to be smoothly distributed, this results in a smooth 
component of convergence given by Kc = ksis — The shear is 7*,c = i^sis = + f^c- In the 
above expressions Da and Dg are the angular diameter distances of the lens and source at redshifts 
of Zd and Zg, and D^s is the angular diameter distance from the lens to the source. Galaxies with 
small velocity dispersions can have > ksis very small radii using this definition, however we 
constrain < ksis for all ^, thus keeping Kc > 0. 



therein) advocates a factor close to unity from dynamical modeling of nearby early-type galaxies and from the 
distribution of observed macrolens image separations. 

^The results are quite insensitive to the star-formation history since most microlensing is due to galaxies at 
redshifts below 1. The exception is in the rate of microlensing of singly imaged quasars by stars in galaxies at higher 
redshifts. 
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Microlensing statistics are described by a combination of the stellar mass-density and the per- 
turbing effects of a continuous component of convergence as well as shear from the galactic mass 
distribution (See Webster et al. 1992 for a detailed description of the dependencies in the case of 
Q2237+0305). Therefore three parameters 7*c) govern the microlensing statistics of any one 

microlensed image. However, to reduce computation, we use a parameter transformation (Paczyn- 
ski 1986; Kayser, Refsdal k. Stabell 1986; Schneider & Weiss 1987) which allows the microlensing 
statistics to be described by an equivalent model having no continuous component of convergence, 
and thus only two microlensing parameters (^,7). Note that this transformation describes a de- 
generacy, and that the dependence due to is merely shifted to other physical parameters. The 
lens equation may be written in component form 

yi,*c ~ (1 l^c l^-k Tc*)2^1,*c 

y2,*c = Kc- K*+ lc*)x2,*c, (4) 

where x = £,/^o,*c and y = rj/rjo^^c are the image and source positions in units of the microlens 
Einstein radius in the lens and source plane. Given point-source magnification /x*c source size 
5*0 the transformations 



K. 



K 



7 



|1 - Kc\ 



7* 



|1 - /^cl 

\2 



/7, = /i*c X (1 - 



S = S^c/^A^-Kc 
rjo = 'no,*c X \/\T^ 

= Co,*c/\/|l - I- 



y = y*c/ \/|i - Kc\ 

X = X*c X - Kc\ (5) 

yield an equivalent lens equation with identical microlensing properties, but no continuous compo- 
nent of convergence, hence 

2/1 = (1 - K- 7)xi 

y2 = {1- K + ^)X2. (6) 

Below we find the differential joint cross-section for k and 7, as well as probabilities for the 

sets of parameters {k, 7, Kc, Zd}. The parameters k, 7 and Hc are not independent quantities 
and the following procedure is used to find First, the probability of an image subject to a 

microlensing optical depth between n and k + Ak at a redshift between and Z(i + /S.Z(i is computed 
(see paper I for details) by taking the derivative with respect to k of the differential cross-section 
for an image being subject to a microlensing optical depth larger than k 

(f'T d dr . , , , 

dndzd dn dz^ 
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The relation ^« = ^{k) is multi-valued, having 2Npairs = 2 or 4 solutions. At each of the 2Npairs 
solutions of = ^{k) we find the constants CV and C2* such that 



Cr (X r? ^ 
dr} 

dr] 



dn 



(8) 



and 

(Cr + C2^) = 1. (9) 

We also find the values of microlensing shear (7^, 72), and the values of smooth matter density 
('^ci' '^02)- ^^is procedure is repeated over a logarithmic grid of values for k, a and z^, with 
spacings Ak, Act and Az^. The resulting sets of values, {n, k^ i, j\, zj} and {k, kJ, 2, 72, -^d} have 
corresponding probabilities 

p({K,7l,4^i,z<i})AKAzd = CV-^^\^=^—AaAKAzd 

and 

pi{K,Y^,i^l^,z})AKAzd = CT^^J^=^—AaAKAzd. (10) 

These probabilities are binned in 7 and integrated over Zd, to find the the differential joint cross- 
section 



l-CO 

dzd / dyH{j,Aj) 
Jo 



dnd-y Jq Jq A7 
where 

^(^•^^' = {o *ele''''^^ <"> 

of observing an image subject to a microlensing optical depth between k and k + An, and mi- 
crolensing shear between 7 and 7 + A7. For impact parameters smaller than the galaxy Einstein 
radius, the SIS lens model produces 2 images. For multiply imaged sources the values of k, 7 
and Kc for the second image are kept in the set in addition to those of the primary image. The 
parameters for the second image are easily computed by noting that the image separation for an 
SIS is always Ax = 2. These additional parameters facilitate calculation of the distribution of total 
magnification for pairs of multiple images (Sec. 3.1). 

The lower rows of Fig. 1 show contour plots of the differential joint cross-section for 
sources that are singly imaged by the galaxy (left), multiply imaged by the galaxy (center), and 
for all images (right). The case shown has a source redshift of = 3 and the contours are spaced 
by factors of \/lO. Several features of the distributions warrant explanation: 
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1) K > 7 for all K, as can be easily seen from their definitions (Eqn. 5). 

2) There is a line of demarcation between regions of k — 7 parameter space accessible to single and 
multiple images. Multiple images formed by an isothermal sphere have ksis = i^* + Kc = 7575 = 
7*c > 1/4. This demarcation line is therefore parameterized by 



+ 3/4I 



and 



Given positive this yields 7 = for k < 1 and 7 < 1/3. 

3) An SIS has ksis = i^-* + i^-c = 'ysis = 7*c = 1/2 at the Einstein radius. The relation between k 
and 7 at the Einstein radius is therefore parameterized by 

and 



\k^ + 1/2\ 

which results in the relation 7 = 1 — k. This is the condition of formally infinite magnification 
required at the Einstein radius. The line l//x = |(1 — k)^ ~ 7^1 = is shown as the dashed line 
in the lower panels of Fig. 1. Note that this line runs through a valley in the probability density, 
which is the result of depletion of images near the galaxies' Einstein radii. The upper panels of 
Fig. 1 show ^ obtained by integrating over 7. The plots clearly show the depletion near the 
Einstein Radius. 



3. The Probability for Quasar Microlensing 

In this section we calculate microlensing statistics by combining probabilities for different com- 
binations of microlensing parameters p{{k, 7, Kc, Z(i})AKAzd with results from numerical microlens- 
ing simulations. A large number (556) of magnification patters were computed over the region of 
interest defined by -^j^ . The computation of these magnification patterns was performed using the 
microlens ray-tracing program, generously provided by Joachim Wambsganss. The values of k and 
7 for which magnification patterns were computed are marked (light dots) in the lower right panel 
of Fig. 1. Each magnification pattern was 25 microlens Einstein radii on a side. One magnification 
pattern of this size does not adequately describe the microlensing statistics for the corresponding 
set of microlensing parameters (Seitz, Wambsganss & Schneider 1994), and as a result probabilities 
calculated for a single set of parameters will not be accurate. On the other hand, our calculations 
of microlensing probabilities average over a large number of patterns having different microlens- 
ing parameters, so that simulation variance will average out. Furthermore, we normalize each 
magnification pattern by the corresponding theoretical magnification {^^h = |(1 — k)^ — 7^|~^). 
Magnification patterns having k < 0.025 were not computed, and statistics of microlensing for 
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(«;, 7) where n < 0.025 were found by multiplying the probabilities computed using the (0.025, 7) 

magnification pattern by (1 — e~'')/(l — e^'^''^^^). We computed 10 realizations of each of the 20 
magnification patterns at k = 0.025, so that statistics for a given shear with k, < 0.025 would not 
be heavily biased by any single pattern. 

In paper I, the Poisson probability that a source will lie inside the Einstein ring of at least one 
microlens was used to approximate the conditional probability distribution for the microlensing 
optical depth near lines of sight to microlensed quasars 

^ = (l_e-)i^. (14) 

anaz dndz 

This approximation has several shortcomings. Firstly, the contribution of 7 to microlensing statis- 
tics is ignored. Secondly, the probability of microlensing at k ~ 1 is assumed to be ~ 1, while 
simulations at k it 7 ~ 1 show that the size of light-curve fluctuations tends to zero (Deguchi & 
Watson 1987; Seitz, Wambsganss k, Schneider 1994). Thirdly, the source size is not considered, 
even though a larger source size results in longer timescales and reduced event amplitudes for mi- 
crolensing events (Wambsganss, Paczynski &: Katz 1989). In addition to these shortcomings, the 
quantity computed by Eqn. 14, namely the fraction of sources microlensed at any one time is not 
necessarily the quantity of interest. Rather it is more useful to know the magnification distribution 
resulting from microlensing, and the likelihood of variability above some threshold level. In the 
following subsections we consider these in turn. 

In what follows, the source quasar is assumed to have a radius of S^^c = lO^^cm (typical scale 
for an accretion disc; found for Q2237-I-0305 by Wambsganss, Paczynski & Schneider 1990 and 
Q0957+561 by Refsdal, Stabell, Pelt & Schild 2000) and to have a uniform top-hat profile. We 
note that the size of the emission region for a fixed observational band will decrease with source 
redshift (for a thermal accretion disc) since wc are observing intrinsically higher frequencies and 
therefore brightness temperatures. Since has been computed for the effective microlensing 
optical depth and shear defined by Eqn. 5, the physical value of source size must be adjusted 
accordingly for each value of Kc- The mass-spectrum is demonstrably unimportant (e.g. Witt, 
Kayser Sz Refsdal 1993), and the microlenses are assumed to have a single mass of O.IM0. 



3.1. Magnification Distributions 

Recall that we find a probability p({k, 7, k^, ^dDA/tAzcj for each set of values {k,'j, k^., z^}. 
To find the magnification distribution due to galaxies containing populations of stars, these prob- 
abilities are convolved with the magnification distribution computed for the corresponding sets 
of microlensing parameters. A coarse magnification distribution ^(k, 7, -s^d) was obtained for 
each parameter set by binning 100 magnifications generated from the magnification map having 
microlensing parameters closest to those of the current set (we also normalize the mean of the 
magnification pattern by the theoretical mean corresponding to the current set of parameters) . For 
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each parameter set the source size was determined in units of effective microlens Einstein radius 

(S/rjo = S^c/'no,*c |1 ~ /^cD- The physical magnifications were then determined from the scaled 
magnifications /i*c = ~ ^c)'^ which were computed by convolving the source profile with the 
magnification map. The resulting distributions were then integrated over Zd and k. 

dr f°° f^" dP 

— =/ dK dZdp{{K,J,Kc,Zd})-^{K,J,Kc,Zd) (15) 

J 10-4 Jo dji 

The procedure misses lines of sight having k smaller than 10~^ (the smallest value considered). 
However wc know that the single image distribution [ ] must be normalized to 1 — Tmult (where 
Tmult is the cross-section for multiple imaging), and that its mean must be (1 — ^Tmuit)/{^ — 'Tmult) 
so that the average of the magnification distribution = ^"^g^" + ^^3^] for all quasars is unity. 
These conditions were fulfilled by adding probability smoothly to the bins between jj, = 0.9 and 
1.1 [wc arc not concerned with the details of the distribution near /x = 1 which must be computed 
from n-body simulations (e.g. Barber, Thomas, Couchman &: Fluke 2000)]. Inaccuracies can 
arise from the integration over the singularity in magnification as a function of n and 7 using a 
finite grid. The multiple image distribution was corrected for this effect by multiplying by the 
analytic SIS distribution over the numerical distribution assuming a smooth mass distribution (i.e. 
no microlensing) . Magnification distributions were computed for = 1, 2, 3 and 4. Example 
distributions for a source at plotted in Fig. 2 for single (^^, left), multiple (^^^^, 

center) and all (^ = '^'^2^^ + ^^^djf^^ right) images. The distribution for multiple images includes 
individual values for both the bright and faint macro-images. Also shown for comparison in the 
central panels is the analytic distribution for a smooth SIS. The figure shows that microlensing 
results in a spread of both the single and the multiple image magnification distributions. Of 
particular note in the multiple image distribution is the small excess of large magnifications over 
the smooth SIS level. In the distribution for single images, we see that the effect of microlensing 
is to create a significant non-zero probability for magnifications with values greater than 2, which 
are not formed by the SIS. 

The distribution (3^^) of total magnifications {utot, the sum of the magnifications for the 
bright and faint images) for multiply imaged quasars was computed using an analogous procedure. 
However instead of computing the distribution of magnifications for a single magnification pattern, 
the distribution was computed for the sum of magnifications computed from pairs of magnifica- 
tion patterns. These magnification patterns had microlensing parameters corresponding to those 
recorded for both images in the sets of parameters described in Sec 2. The resulting distribution 
is shown in Fig. 3 for a source at Zs = 3. The analytic distribution for a smooth SIS is again 
shown for comparison. We find an excess of large magnifications and a non-zero probability for 
lifot < 2, the minimum value for the total magnification of a smooth SIS. In Sec. 3.3 we will use this 
distribution to compute the effect of microlensing on the magnification bias for multiple imaging 
in optical gravitational lens surveys. 



- 11 - 



3.2. The Probability of Quasar Microlensing Variability 

In this subsection we compute the probabihty that the magnitude of a quasar wih vary by 
more than Am due to microlensing during a 10 year period. For this calculation we assume that 
the velocity components v^, and Vy of the galaxy {vd), source {vg) and observer [vo) are Gaussian 
distributed with cr^ = 4Q{)kmsec~^. From Kayser, Refsdal k, Stabell (1986) we compute the 
effective source plane transverse velocity 



Vpff*r = \lv'^tf +v'^tt , where 

Veff,*c,x = T-, FT + T"; TT + T"; 

1 + ZdDd 1 + ZdDd l + Zg 

= TT7d^ + TT7dD-d + TT^- ^''^ 

Since the probabilities p{{k,j, Kc, Zd})AKAzd have been computed in terms of the effective mi- 
crolensing optical depth and shear defined by Eqn. 5, the physical value of transverse velocity must 
be adjusted accordingly (Eqn. 5) for each value of Kc 



'^eff = Veff,*c 



/VlT^. (17) 



3.2.1. microlensing variability cross-sections 

Again recall that we find a probability p{{k, 7, Kc, Zd})AKAzd for each set of values {k, 7, Kc, Zd}- 
To find the conditional joint differential cross-section for values of k and 7 near lines of sight to 
microlensed images, these probabilities were multiplied by the probability / that the source will 
be microlensed. In this case we take f{Am\K,^,Kc,Zd) to be the fraction of light-curves m{t) 
with 10 year monitoring periods that vary by more than Am (note, in paper I / was taken to 
be 1 — e"'^, Eqn. 14). We computed f{Am\K,j,Kc,Zd) for each parameter set from 100 light- 
curves generated using the magnification map having microlensing parameters closest to those of 
the current set (we also normalize the mean of the magnification pattern to that of the current 
set of parameters). The light-curves had directions perpendicular, parallel and at 45 degrees to 
the shear with weightings of 0.25, 0.25 and 0.5 respectively which approximates random directions 
for the transverse velocity. For each parameter set the transverse velocity in units of Einstein 
radii per second is = -^^^ji^^- The true magnification is related to the magnification 
computed from the integral over source size by ij,^f. = /x/(l — Kc)^, and the resulting light-curve is 
m{t) = — 2.51 log(/x*c) + const. Values of f{Am\K,j,Kc,Zd) x p{{k,^, Kc, Zd}) were then binned in 
7 and integrated over Zd as before to yield the conditional differential joint cross-section for values 
of K and 7 near lines of sight to microlensed images 

d^TML r" , , frrr A x/(^"^l'^>7,«c,^d) X p({K,7',Kc,2d}) 



r dzd r dl'Hi^, ^^^f{Am\K,j,K,,Zd)xp{{K,j,n,,Zd})_ ^^^^ 
Jo Jo ^7 
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The lower panels of Fig. 4 show for the case of a source at = 3 and Am = 0.5 for singly 

imaged (left panel), multiply imaged (considering each image separately, center panel) and all 
(right panel) sources. The distribution shows many of the features seen in Fig. 1, namely depletion 
near the Einstein Radius, demarcation between regions having singly and multiply (by the galaxy) 
imaged sources (light dashed line), and k < 7 for all k. However, the distribution is suppressed in 
the regions of both low and high k. The upper panels of Fig. 4 show obtained by integrating 
over 7. Distributions are plotted for Am = 0.5, 1.0 and 1.5 magnitudes. The integrals of 
these distributions TML(Am,Zs), which we term the microlensing cross-section, i.e. the fraction of 
sources that vary by more than Am per decade due to microlensing are also given. 

The microlensing cross-section rML(Am, Zs) is plotted as a function of Am in Fig. 5 for Zg = 1, 
2, 3 and 4. The left, central and right panels show values for single quasar images, multiple quasar 
images (considering both images separately) and all images respectively. For example, ~ 1 in 1000 
quasar images at = 3 vary by more than Am = 0.5 magnitudes due to microlensing. the majority 
of this microlensing occurs in sources that are also multiply imaged by the galaxy. As expected, 
large amplitude variability is rarer than low amplitude variability, and there is a 10 fold decrease 
in rate between Am = 0.5 and Am = 2.5 magnitudes in all examples. Microlensed variability at 
all levels is more likely at higher source redshift. 

Microlensing cross-sections were also computed under the assumption that the dark matter is 
entirely composed of compact objects. The results are shown as the grey lines in Fig. 5. The SIS does 
not have finite mass, and so the surface mass distributions were truncated at an outer radius ^max 
(proportional to for each galaxy) such that the sum of the masses constitutes the cosmological 
density CI. The numerical procedure was checked by noting that the mean of the distribution 
should equal the result obtained for the optical depth of randomly distributed objects (Press & 
Gunn 1973; Turner, Ostriker &; Gott 1984). Since much of the mass inside the galaxies Einstein 
ring is in stars, the microlensing cross-section in multiple images is quite insensitive to the addition 
of dark compact objects, giving an enhancement of only a factor of ~ 2 — 3. On the other hand, 
while the microlensing cross-section due to stellar populations is very small for single images due 
to the rapid decline in the density of stars beyond the critical radius for multiple imaging, most of 
the dark matter mass lies beyond this radius. As a result, the inclusion of dark compact objects 
boosts the single image microlensing cross-section significantly to ~ 10%. This is expected since 
in the formalism of Press & Gunn (1973) we find that the optical depth is a few tenths of Q. We 
note that this is a very coarse approximation to lensing by the large scale dark-matter distribution, 
which is a topic beyond the scope of this paper. However, an unavoidable consequence of Fig. 5 
is that comparison between variability of multiply imaged quasars, and quasars with lines of sight 
near galaxies will be a powerful and achievable method for detection or rejection of compact objects 
as dark matter candidates in the halos of galaxies. 

The microlensing cross-sections described above assume a 1-d dispersion in the proper motions 
of galaxies of i; = 400fcmsec~^, microlens masses of (m) = O.IMq and a source size of S* = lO^^cm. 
However we can discuss the results qualitatively for other choices for these parameters by noting 



-13- 



that event peak amplitudes scale with (m)^/^ and S"^/^, while event rates scale with (m)~^/^ and 
V (Witt, Kayser & Refsdal 1993). Using these relations we approximate the general results for the 
microlensing cross-sections shown in Figs. 5 and elsewhere in the paper for a 1-d proper motion 
velocity dispersion of v' , source size S' and microlens mass {m)' as 

where 

Am ~ Am' — 2.51 log 

As an example, for a source S' = lOS microlensed by microlenses of mass (m)' = O.OIMq wc 
find that T']^j^{Am') 3.16rML(Am = Am' + 1.88). For Am' = 0.5 magnitudes, this results in 
T'j^]^{Am') ~ 4 X 10""^ or half the value of TML{^m = 0.5). 



(m)' 
O.IMq 



S' 



1015 



cm 



(19) 



3.2.2. the probability of microlens redshift 

The probability distribution for the redshift of galaxies whose stellar populations result in 
the microlensing of back-ground sources was computed using an analogous approach. The resulting 
distributions ^^g^ for the case of = 3 are plotted in Fig. 6 for singly imaged (left panel) , multiply 
imaged (center panel) and all (right panel) sources. Distributions are plotted for Am = 0.5, 1.0, 
1.5, and 2.5 magnitudes, and can be compared to the distribution of macrolens redshifts (light 
line). The distribution describes a low and narrow range of redshifts of microlensing galaxies that 
also produce multiple imaging. This peak is bounded from below by the low number of macro-lens 
galaxies at low redshift (e.g. Turner, Ostriker & Gott 1984). From above the peak is bounded 
by two factors of Dg/D^, one for the size of the source with respect to the microlens Einstein 
radius projected into the source plane (i.e. larger fluctuations for larger Ds/D^), and one for the 
timescale which decreases in proportion to Dg/Dd (resulting in an increase in microlensing rate). 
Furthermore, at higher redshifts the typical impact parameter for multiple images is larger with 
respect to the scale radius Rq of the stellar distribution. Multiple images are therefore typically 
subject to lower values of k^. As a result of these factors, the typical redshift of microlensing galaxies 
is lower than the typical macrolens galaxy redshift. While the microlensing cross-section at higher 
lens galaxy redshifts for singly imaged quasars is decreased due to the factors of Dg/D^ already 
mentioned as well as from the reduced microlens population above ~ 1, the lower macrolensing 
cross-section for high-redshift lenses gives more single images lines of sight through regions of high 
Ki, in those galaxies. As a result is broad for single images, in contrast to results for multiple 
images. Note however that the larger amplitude variability in single images is still only found for 
galaxies at low redshift since large amplitudes require a source that is small with respect to the 
caustic structure. 
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3.2.3. observable microlensing statistics 

The results thus far have referred to probabihties for microlensing at fixed source red-shift, 
and have not included the effects of magnification bias. However, both the distribution of source 
redshifts, and the magnification bias need to be considered for a connection to be made with ob- 
servations. In this subsection we find observed microlensing probabilities by combining differential 
cross-sections with the empirical double power-law quasar luminosity function Zs) of Pci (1995) 
(based on data presented by Hartwick &: Schade 1990 and Warren, Hewctt &; Osmer 1994) and 
the break luminosity evolution described by Madau, Haardt & Rees (1999). To convert the em- 
pirical luminosity function back to relative number counts above a limiting apparent B-magnitude 
(m^) we convert ms to a limiting intrinsic luminosity using the luminosity distance and 
a /c-correction found from the procedure described in Pei (1995) and M0ller & Jakobsen (1991). 
The magnification due to gravitational lensing allows observations at a fixed limiting magnitude 
to reach further down the luminosity function and include more sources (e.g. Turner 1980). The 
resulting bias for the fraction of sources exhibiting the lensing phenomenon of interest often plays 
an important role in the analysis of lens samples. Magnification bias for a varying source requires 
careful treatment that depends in detail on the observed sample. In the remainder of this subsec- 
tion we describe two different examples of the statistics of microlensing variability in hypothetical 
monitoring campaigns. 

In the first example we compute the fraction of quasar images microlcnscd by more than Am 
per decade in a continuous blind survey to a fixed limiting magnitude tub- This type of survey will 
be performed as part of the Sloan Digital Sky Survey (York et al. 2000) where a stripe of sky will be 
repeatedly scanned. Let Psgie{zs) / Pmulti^s) be the fraction of quasar images that are microlensed 
and in single / multiple image systems (at source redshift Zg). These are given by 

P{zs, Am, mum) = , , ^ . (20) 

where Pl^t{zs, Am,miira) = P^gi^^Zs, Am,mum) + P^^i^{zs, Am^mum)- The quantities P'{zs) = 
t{zs)B{zs) (where B is the magnification bias) are the fraction of images that are are in single 
or multiple images systems (images in multiple image systems are considered separately). The 
probability that an image is detected during the survey is the cross-section t{zs) multiplied by the 
bias calculated at the maximum light-curve magnification fXmax- Since this maximum is subject to 
microlensing, the bias is given by the expression 



B{zs,Am,miim) = 
1 



t{Zs) 



oc rZs 

dK / dZdp{{K,J,Kc,Zd}\Zs 
^0 



dHmaxd^{t^^l^Zd,Zs)N{> ^,Zs) 
N{> Lum,Zs) 



(21) 



where N(> Lum, Zg) = dL^iL, Zs). The detection of variability larger than Am requires that 
the image be detectable when Am magnitudes fainter than the light-curve maximum. Therefore, 
while the bias for the detection of an image is computed using the magnification at the light-curve 
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maximum, the bias for detection of a microlensed image should be calculated using a magnification 
of 

A*i)ar — Mmaa; 10 The parameter F in Eqn. 20 is the fraction of images in single / multiple 

image systems that are microlensed by more than Am. This is given by 



F{zs,Am,mHmj 

— — / dn I dZdP{{n,-i,Kc,Zd}\zs) 
yZs) Jo Jo 



^(Zs) 



dP , f{iJ.max\^rn)N{> 



p dP 

Jo ^l^max 



f-^var ' 



^ fJ'max ' 



(22) 



where /(//mail^^) is the fraction of light-cTirvcs with [imax that show variability above the level 
Am. The integrals over iimax in Eqns. 21 and 22 were computed via Monte-Carlo. Integrating 
Eqn. 20 over source redshift we find 



/o°° dzs ^N{> Lu,n, Zs)P{zs, Am , m^m) 

I^dZs'^N{> Lum,Zs) 



P(Am, mum) = (23) 



where Vcm is comoving volume. For this calculation P{zs) was computed for Zg = 1, 2, 3 and 

4 (by Zs = 4: Lya has moved through B-band), and interpolated in Zg during convolution with 
N{> Liira,Zs). The resulting values of Psgie{Am,mii,n) (left), Pmulti^fn^mu^) (center) and 
Paii{Am,miirn) = Psgie{Am, mum) + Pmulti^m, miim) (right) arc plotted as a function of Am 
for several values of m^m in Fig- 7. Also shown in Fig. 7 (grey lines) are the corresponding results 
for halos that are also composed of dark compact objects. 

Magnification bias introduces two competing effects with regard to the microlensing rate. 
First, low magnification images near the center of the lens, and those outside the Einstein radius 
are less likely to be observed, while the probability of observing images near the Einstein radius 
having values of k favorable to microlensing variability is enhanced. This leads to an increase in 
the microlensing rate. However magnification bias also reduces the fraction of images that are 
observed to undergo large fluctuations. The bias described in Eqn. 21 suppresses the frequency of 
observed large amplitude fluctuations because of the larger difference between the number of sources 
detectable during the survey at light-curve maximum and the number detectable while Am below 
the light-curve maximum. The net result is to increase the rate of stellar microlensing variability at 
a level Am > 0.5 by up to a factor of 10 (for mu^n = 17) in images of both macro-lcnscd and singly 
imaged quasars. On the other hand, when we consider magnification bias in simulations where dark 
matter is composed of compact objects, wc find the microlensing rates are not increased relative to 
the no bias case. In the case of multiple images, moderate optical depth is present right out to the 
critical radius for multiple imaging, so that there is no gain in terms of a more favorable k. Thus 
the microlensing rate decreases with increased magnification bias, which is the opposite behavior 
to that of stars. 



For our second example, we compute the rate of microlensing in the images of lensed quasars 
that have been previously identified in a survey for multiply imaged gravitationally lensed quasars. 
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In this case the bias corresponding to the microlensing rate for images formed at a given set of 
microlensing parameters (k, 7) is simply the bias for observing an image with that set of parameters 
in the initial survey. In general, this bias is computed using the sum of the magnifications of multiply 
imaged sources {/J^tot) since the resolution of survey quality data is lower than that required to 
identify individual lensed images (e.g. Webster, Hewitt Sz Irwin 1988). We compute the fraction of 
quasar images in macrolensed systems that exhibit microlensing. This is the most readily observed 
statistic since (following determination of a time delay) microlensing variability can be distinguished 
from intrinsic variability which is seen in all images. The fraction (.Fmw/t) of macro-images (of 
quasars at redshift Zg) to undergo microlensing variability is 

Fmult{zs,^m,mum) = 

i^c,Zd}\zs) diJ,tot-;S^{i^,l,Zd,Zs)fiAm\K,-f,Zd,Zs)N{> ^,Zs) 

-(24) 





[Jo dntotdftoSi^^^^^d,Zs)f{/^m\K,-i,Zd,Zs)N{> ^';;^,^s)J 


dK J^' dZdp{{K, 7, Kc, Zd}\Zs) 


' d^^tot^S^,^,Zd,Zs)N{> ^,zs)_ 





where /(Am|K;, 7, z^, z^) is the fraction of light-curves that vary by more than Am magnitudes, 
and -^^{f^il-, Zdi Zg) is the normalized probability distribution for the sum of image magnifications 
(Sec. 3.1). Integrating Eqn. 24 over source redshift we find 

^ ,^ , dZs^^N{> Liirn,Zs)Fmult{Zs,Am,miim) ^^^^ 

Values of Fmuit are plotted as a function of Am for several values of mum in Fig. 8. Also shown 
in Fig. 8 (grey lines) are the corresponding results for halos that are composed of dark compact 
objects. 

We find that 30 — 50% of multiply imaged quasars should vary by more than Am = 0.5 
magnitudes during a 10 year period. This rate drops by a factor of 5 for large amplitude fiuctuations 
which is shallower than the previous example since the bias is not calculated from the light-curve. 
Note that these rates imply that microlensing induced variability is not uncommon in macrolensed 
quasars, even though the Einstein radius crossing time is many decades for most redshifts. This is 
because in the high k and 7 environments where most macro-images are found the caustic network 
has a typical scale- length significantly less than Furthermore, only part of a source need cross 
a caustic for significant variability to be observed. While the rate for microlensing by stars is 
very sensitive to magnification bias, the rates in the presence of dark compact objects in the dark 
matter halo are quite insensitive to the limiting magnitude. A large magnification bias results in 
an increased fraction of images near the galaxies Einstein radius where is typically around a 
10th. Thus the microlensing rate is increased. On the other hand if the dark matter is composed 
of compact objects, the microlens optical depth is always greater than 0.25 for multiple images. 
Hence magnification bias does not make much difference in this case. The rate does not drop 
appreciably with magnification bias, since unlike the previous example the bias was calculated 
independently of the light curve. The situation is described by Fig. 9 which shows the observed 
probability distribution for k assuming different limiting magnitudes. The left hand panel of Fig. 9 
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shows distributions assuming microlensing by stars, while the right hand panel assumes the halo 
is also composed of compact objects. Magnification bias imposes a sharp peak on the distribution 
near k = 0.1 for stars, increasing the microlensing rate since multiple images are less likely to be 
observed at low k. 



3.3. The effect of microlensing on magnification bias and the distribution of fiux 

ratios 

We have seen that microlensing by stars affects the magnification distribution for a quasar 
image, as well as the distribution for the sum of the magnifications of multiply imaged quasars. 
Since microlensing in the two images is independent, we also expect microlensing to affect the flux 
ratios of images in macro-lensed systems. Furthermore, if microlensing affects the magnification 
distribution, then it might also qualitatively affect calculations of magnification bias for the fraction 
of multiply imaged quasars. As mentioned in the Introduction, Bartelmann &; Schneider (1990) 
computed the effect of microlensing on the magnification bias and macro-image flux ratios. Under 
the assumptions of a point source, analytic forms for the magnification distributions (computed 
for different values of k and 7), and SIS galaxies having a fraction of density in compact objects 
with an outer radius, they concluded that microlensing could severely affect both the magnification 
bias and the flux ratio distribution. We are now in a position do better on all these scores, having 
computed magnification distributions numerically, for finite source size, and for a population of 
microlenses that have distributions resembling known stellar populations. 

Using the luminosity function of Pei (1995) we compute the magnification bias at source 
redshift Zg 

Biz ) = - ' ^ - (26) 

Jo a^i [Tmuityzg) ^^^^^ \^l.tot=^i' ^ d/x' i^^y^ hm'^' 

which results in an overall bias of 

^ ^ dZsT{Zs)^N{> Lu„.,Zs)B{Zs) 

Note that Eqn. 26 includes the possibility that singly imaged quasars are magnified by microlenses 
beyond the galaxy Einstein radius. Fig. 10 shows the magnification bias for multiple imaging as a 
function of the limiting magnitude. The bias for a smooth SIS is shown for comparison (dot-dashed 
line). There is a slight increase in the bias if microlensing is assumed, however the effect is very 
small. Apparently additional bias due to the small increased probability for high magnifications is 
balanced by the possibility of having ^tot < 2, as well as by the increased number of singly imaged 
sources (Bartelmann &; Schneider 1990). At redshifts where macrolensing is most likely, the distance 
ratio is of order 1, and so the source size is comparable to the projected size of the microlens Einstein 
radius. In this case the magnification distribution is quite narrow (Wambsganss 1992). Hence the 
excess of high magnifications formed by microlensed images is not as large as those predicted by the 
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point source distribution, and as a result the effect on the macrolensing magnification bias is small, 
in contrast to the findings of Bartclmann &: Schneider (1990). A similar result is found assuming 
dark matter to be composed of compact objects. 

Finally, we compute the distribution of flux ratios R = ^ (where is now defined to be the 
image with the larger impact parameter, i.e. 1 < xi < 2 or 1/2 > ki > 1/4). The distribution 
is shown in Fig. 11 for a source at Zg = 3. The analytic distribution for a smooth SIS is shown 
for comparison. As a result of microlensing the likelihood of the distribution mode is lowered by 
a factor of ~ 2, and a small excess of large flux-ratios is formed. This quantifies the often made 
statement that microlensing will influence flux ratios, and illustrates why optical flux ratios should 
not generally be used as constraints for models of gravitationally lensed galaxies. 



4. The Probability of Microlensing for Gamma Ray Burst Afterglows 

Just as microlensing promises to probe the central engines of quasars, microlensing of GRB 
afterglows offers a means to probe their structure (Loeb & Perna 1998). Recently an anomalous 
event of 0.95 magnitudes was observed in the afterglow light-curve of GRB000301C and interpreted 
by Garnavich, Loeb and Stanek (2000) as showing features consistent with the microlensing hy- 
pothesis (see also Gaudi, Granot & Loeb 2001). With the afterglow light-curve of GRB000301C 
as motivation Koopmans &: Wambsganss (2000) determined the a-posteriori probability that the 
afterglow was microlensed and found the probability to be small but not prohibitive. Furthermore, 
Mao & Loeb (2001) demonstrate that microlensing at the few percent level should be observed 
in all afterglows on timescales of a year (although they neglect clustering of microlenses) . In this 
section we extend the calculation of Koopmans and Wambsganss (2000) and consider the stellar 
population as microlenses (both alone and in addition to dark compact objects in the halo), as well 
as the fraction of microlensed GRB afterglows that will also be multiply imaged. 

Our calculation follows the method described in Sec. 3.2.1. However the GRB afterglow source 
is assumed stationary, while the relative motion is provided by its expansion. We assume the 
afterglow to be described by an expanding ring of constant surface brightness. The ring has a 
radius R{t) = Ro{t/days)^^^ that expands as a function of time, and a width oiW x R where i?o is 
the width on day 1 (Waxman 1997). We take i?o = 3 x 10^^(1 -|- Zs)~^cui and W = 0.16 (Loeb & 
Perna 1998). The ring radius R{t) must be adjusted for the smooth matter correction (Eqn. 5) as 
before. We assume that the afterglow is monitored for 30 days (in the observers frame) following 
the burst. The luminosity function of GRB's from which the afterglows arc identified is poorly 
known, and we do not consider magnification bias in the calculations of microlensing variability 
cross-sections. 

The lower panels of Fig. 12 show '^J^^ , the conditional joint cross-section for values of n and 
7 near lines of sight to microlensed images of GRB afterglows for the case of a source at Zg = 3, 
and Am = 0.5 for singly imaged (left panel), multiply imaged (center panel, images are considered 
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separately) and all (right panel) sources. Note that values of for k < 0.025 were calculated 

by multiplying ^^^^|«;=o.025 by (1 — e~'^)'^/{l — e"^-'^^^)^ (where the square is included because 
the motion is in both dimensions). The distribution shows many of the features seen in Figs. 1 
& 4, including depletion near the Einstein radius, demarcation between regions having singly and 
multiply (by the galaxy) imaged sources, and k < 7 for all k. The upper panels of Fig. 12 show 
obtained by integrating ^^^^ over 7. Distributions are plotted for Am = 0.5, 1.0, 1.5, and 
2.5 magnitudes. The microlensing cross-sections rML(Am, z^) are also given. 

We found probability distributions for the redshift of galaxies whose stellar populations result 
in the microlensing of back-ground GRB afterglows. The resulting distributions for the case 
of = 3 are plotted in Fig. 13 for singly imaged (left panel), multiply imaged (center panel) 
and all (right panel) sources. As before distributions are plotted for Am = 0.5, 1.0, 1.5, and 2.5 

magnitudes, and the distribution of macrolens redshifts is included for comparison (light line). The 
shape of these distributions are similar to those for quasars. However, the distribution for singly 
imaged GRB afterglows is more peaked than that for quasars, and has a faster decline with redshift. 

We have computed the microlensing cross-section TML{Am, Zg) for GRB afterglows as a func- 
tion of Am and plotted the results in Fig. 14 for Zg = 1, 2, 3 and 4. The left, central and right 
panels show values for single, multiple (by the galaxy) and all images respectively. There are several 
points of interest. Firstly, the relative rates of microlensing between singly imaged and multiply 
imaged sources are similar for GRB afterglows and quasars. Around 1 in 1000 Zg = 2> afterglows will 
be microlensed (Am = 0.5 magnitudes) by stars, however nearly all of these will be singly imaged. 
Secondly, the microlensing rate falls off steeply with the amplitude of the fluctuation Am, dropping 
by a factor of ~ 20 between Am = 0.5 and Am = 1.5. This is steeper than the corresponding 
dependency for quasars (see Fig. 5) . The largest fluctuations require a source that is small with re- 
spect to the projected Einstein radius, i.e. a low redshift for the galaxy containing the microlenses. 
For quasars this is accompanied by an increase in the effective transverse velocity, which helps to 
balance out the reduced probability for intercepting a caustic. However, the expansion velocity of 
the GRB afterglow is independent of z^, hence the more rapid decline of microlensing rate with 
Am. Fig. 14 also shows results obtained under the assumption that dark matter is composed of 
compact objects (light lines). We see an increase of a factor of ~ 10 for macro-lensed afterglows, a 
larger discrepancy than for quasars. However, the fraction for all images (AM > 0.5 magnitudes) 
rises to more than 1 in 10. Note that we get a rate of ~ 1 — 3% {zg > 2) for Am > 0.95 magnitudes 
which is consistent with the value of 5% found by Koopmans & Wamsganss (2001) for an optical 
depth of ~0.25 over the whole sky. Features of the sort seen in the afterglow of GRB000301C 
should therefore be common if dark matter is composed of compact objects. This suggests that the 
large numbers of afterglows that are expected to be discovered by the upcoming swift satellite will 
make an ideal probe of the cosmological density in dark compact objects since even a moderate 
number of featureless light-curves will provide a tight constraint. Furthermore, the high probability 
for GRB microlensing allows for the possibility of probing the distribution of compact objects in 
the halos around galaxies. 
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Finally we determine the fraction of multiply imaged GRB afterglows that are subject to 
microlensing. The results are shown in Fig. 15. Only around 1 in 10 macrolensed GRB afterglow 
images will vary by more than 0.5 magnitudes due to stellar microlensing. This rate rises to around 
60% if the dark halos of galaxies are also composed of compact objects. 

5. Conclusions 

We have computed cosmological microlensing statistics for the populations of stars in ellipti- 
cal/SO galaxies. The stars were distributed with de-Vaucoleurs' profiles and embedded in a smooth 
dark halo to form an overall singular isothermal density distribution (SIS). Our calculation com- 
bines the joint cross-section for the microlensing optical depth k and shear 7 with a large number of 
numerical microlensing magnification patterns, which we assume resulted from microlenses each of 
O.IMq. Wc calculated the fraction of quasar (with assumed source size of lO^'^^cm) images that vary 
by more than Am magnitudes during 10 years of monitoring. We find that most microlensing will 
be observed in multiple images, and that the few cases of microlensing variability in sources that 
are not multiply imaged will primarily have 7 > 0.1 where the point-mass lens does not provide an 
adequate description. Furthermore, the majority of microlensed multiple images of quasars will be 
due to a lensing galaxy at a redshift lower than expected for macrolensing. In contrast, the rate 
of microlensing for sources that arc not multiply imaged by the intervening galaxy is reasonably 
insensitive to the lens-galaxy redshift. 

Wc computed microlensing variability rates for two different examples of hypothetical surveys, 
first, for quasars in a continuous survey of a region of sky to a fixed limiting magnitude. For a 
limiting B-magnitude of iub = 21 we find that 1 quasar in 500 should vary by more than 0.5 
magnitudes during 10 years of monitoring due to microlensing by stars. However 90% of these are 
in multiply imaged systems. Macrolensed quasars therefore dominate microlensing statistics. If the 
dark halo (truncated so that the total mass density equals the critical density) is also composed of 
compact objects, then the fraction of quasar images subject to microlensing variability larger than 
0.5 magnitudes rises to ~10%. On the other hand, the number of multiply imaged microlensed 
sources is quite insensitive to the inclusion of the additional compact objects. Therefore, if dark 
matter is composed of compact objects (as opposed to stars supplying all the microlenses) then 
microlensing is 100 times as common and is dominated by singly imaged quasars. The comparison 
of variability rates of lenscd and unlensed quasars will therefore provide a powerful probe of the 
existence of dark compact objects. For our second example, we assumed that a sample of multiply 
imaged quasars had been previously selected as having the sum of the macro-images brighter than 

= 21. In this case we find that 1 image in 3 multiply imaged quasars should vary by more than 
0.5 magnitudes during 10 years of monitoring. 

We have computed magnification distributions for individual quasar images in both single and 
multiple image systems, the magnification distribution for the sum of macrolensed images and the 
distribution of flux ratios for multiple image systems. Microlensing results in a spreading of all these 
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distributions. In particular, in the presence of microlensing single images can have magnifications 
greater that 2, and the bright image of a multiply imaged quasar can have a magnification smaller 
than 2 (neither is formed by the SIS). From the distribution for the sum of the image magnifications 
we computed the magnification bias for the discovery of multiply imaged quasars that results 
from the inclusion of microlensing in the magnification distribution. Surprisingly, the effect on 
magnification bias for lens surveys is very small. However, we find that the inclusion of microlensing 
has a significant impact on the distribution of flux ratios, resulting in an excess of large values, and 
lowering the likeli-hood of the mode by a factor of 2 with respect to the smooth SIS distribution. 
This illustrates the point that optical flux-ratios make poor constraints for lens models. 

We have also computed the fraction of GRB afterglow light-curves that exhibit microlensing. 
Specifically, we calculated the fraction of GRB afterglows that are perturbed by more than Am 
during the first 30 days. Qualitatively the results differ from those for quasars in that large ampli- 
tude fluctuations are much less common with respect to small amplitude fluctuations. We find that 
only 1 GRB afterglow in 1000 will vary by more than 0.5 magnitudes due to microlensing by stars. 
However most of these will also be multiply imaged by the galaxy. If the halo is comprised of dark 
compact objects we find that the fraction of GRB afterglows that vary by more than AM > 0.5 
magnitudes rises to more than 1 in 10. Features of the sort seen in the afterglow of GRB000301C 
(but with a smaller amplitude) should therefore be very common if dark matter is composed of 
compact objects. If a GRB afterglow is multiply imaged, then we find that there is ~1 chance in 
10 that microlensing by stars will also be observed at a level greater than 0.5 magnitudes, but a 
60% chance if the dark matter is also in compact objects. 
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Fig. 1. — Lower Rows: Contour plots (with contours spaced by multiples of \/lO) for the joint 
differential cross-section of k, and 7. The dashed lines show the condition for infinite magnification 

1 — K = ±7, and the dots in the right panel the parameters of the computed magnification maps. 
Upper panels: Differential probability for k. The left, center and right panels show probabilities 
for single, multiple and all images respectively, the source redshift was Zg = 3. 
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Fig. 2. — Magnification distributions for quasar images. The left, center and right panels show 
distributions for single, multiple and all images respectively. The dot-dashed lines in the central 
panel show the distribution for a smooth isothermal sphere and the source redshift was Zg = S. 




Fig. 3. — Magnification distributions for the sum of images in multiply imaged quasars. The 
dot-dashed lines show the distribution for a smooth isothermal sphere and the source redshift was 
Zs = 3. 
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Fig. 4. — Lower Rows: Contour plots (with contours spaced by multiples of VTO) for the joint 
differential cross-section for k and 7 near the lines of sight to microlensed (Am > 0.5) quasars. 
The dark dashed lines show the condition for infinite magnification 1 — k = ±7, and the light 
dashed lines separate regions of single and multiple images. Upper panels: Differential probability 
for K near the lines of sight to microlensed quasars. The solid, dashed, dot-dashed and dotted lines 
correspond to Am > 0.5, 0.75, 1.5 and 2.5. The left, center and right panels show probabilities for 
single, multiple and all images respectively, the source redshift was Zs = 3 



-29- 




Fig. 5. — Plots of the quasar microlensing cross-section tml verses AM. The left, center and right 
panels show values for single, multiple, and all images respectively. The dark lines correspond to 
microlensing by stars, while the light lines assume dark matter to be in the form of compact objects. 
Values are shown for = 1, 2, 3 and 4. 
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Fig. 6. — Probability distributions for the redshift of galaxies responsible for microlensing of 
quasars. The left, center and right panels show values for single, multiple and all images re- 
spectively, and the solid, dashed, dot-dashed and dotted lines correspond to Am > 0.5, 0.75, 1.5 
and 2.5. The light line in the central panel shows the distribution of macrolens redshifts. The 
source redshift was Zg = S. 
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Fig. 7. — Plots of the probability of microlensing Pml verses AM in a blind monitoring campaign. 
The left, center and right panels show values for single, multiple and all images respectively. The 
dark lines correspond to microlensing by stars, while the light lines assume dark matter to be in 
the form of compact objects. Values are shown for limiting B-magnitudes of = 17, 19, 21 and 
23. 




Fig. 8. — Plots of the probability that a macrolensed image will be microlensed {F^ult) verses AM 
for monitoring of known lenses. The dark lines correspond to microlensing by stars, while the light 
lines assume dark matter to be in the form of compact objects. Values are shown for limiting 
B-magnitudes of m^m = 17, 19, 21 and 23. 
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Fig. 9. — Differential probability of k for multiply imaged sources in the presence of magnification 
bias. The left hand panel corresponds to microlensing by stars, while the right hand panel assumes 
dark matter to be in the form of compact objects. Distributions are shown for limiting B-magnitudes 
of miim = 17, 19, 21 and 23. 




Fig. 10. — Plots of the magnification bias B verses limiting B-magnitude mum The solid and 
dot-dashed lines correspond to biases including and not including microlensing. 
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Fig. 11. — Flux ratio {^1/^2) distributions for the images in multiply imaged quasars. The dot- 
dashed lines show the distribution for a smooth isothermal sphere and the source redshift was 
Zs = 3. 
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Fig. 12. — Lower Rows: Contour plots (with contours spaced by multiples of \/lO) for the joint 
differential cross-section for k and 7 near the lines of sight to microlensed (Am > 0.5) GRB 
afterglows. The dashed lines show the condition for infinite magnification 1 — k = ±7. Upper 
panels: Differential probability for k near the lines of sight to microlensed GRB afterglows. The 
solid, dashed and dot-dashed lines correspond to Am > 0.5, 0.75 and 1.5. The left, center and 
right panels show probabilities for single, multiple and all images respectively, the source redshift 
was Zs = 3. 
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Fig. 13. — Probability distributions for the redshift of galaxies responsible for microlensing of 
GRB afterglows. The left, center and right panels show values for single, multiple and all images 
respectively, and the solid, dashed and dot-dashed lines correspond to Am > 0.5, 0.75 and 1.5. The 
light line in the central panel shows the distribution of macrolens redshifts. The source redshift 
was Zs = 3 . 




Fig. 14. — Plots of the microlensing cross-section tml verses AM. The left, center and right 
panels show values for single, multiple and all images respectively. The dark lines correspond to 
microlensing by stars, while the light lines assume dark matter to be in the form of compact objects. 
Values are shown for = 1, 2, 3 and 4. 
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Fig. 15. — Plots of the microlensing cross-section divided by the muhiple imaging cross-section 
(TML/Tmult) vcrscs AAI . The dark hnes correspond to microlensing by stars, while the light lines 
assume dark matter to be in the form of compact objects. Values are shown for Zg = 1, 2, 3 and 4. 



